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(3) The Net, trained on the 72 realisations as in (1), was then shown the single 'true' IRAS 
observation. To our surprise the Net assigned 100 % probability of belonging to a CDM model (and 
% for the other two models). It was puzzling that the Net favoured one model so strongly. But 
then we realized that one of the 72 realisations happened to have Harmonic coefficients very similar 
to those of IRAS. On one hand it indicates the strength of the Net in recognizing similar objects. But 
it raises the question of 'fair statistics', i.e. how common such a simulation is. 

(4) To cope with the problem of the small number of simulations we generated 10 Bootstrap 
combinations of realisations (allowing repetition). We trained the Net on each of them and presented 
each trained Net with the IRAS Harmonics. In some cases the Net favoured a CDM model, in other 
cases it favoured n = — 1 model, but rarely the Poisson model. The average (with large scatter) 
is 27% 45% and 17% for CDM, n = — 1 and the Poisson models, respectively. The fact that the 
probabilities do not sum up to 100% indicates that our Net is not yet perfect, due to the small number 
of simulations. 

These results, hinting at excess power relative to standard CDM, are in accord with more conven- 
tional statistical method of Maximum Likelihood for the Harmonic amplitudes ([24], [25]) and other 
studies (e.g. [5]) To explore to what extent the Net actually used the phase information we also calcu- 
lated the cross-correlation (a^iRAS a Tsim) an< ^ f° un( i, not too surprisingly, that for / < 4 (large scales) 
they are consistent with random phases (which were the initial conditions of the simulations). Hence 
in this example it is likely that the Net mainly made inference on the basis of the amplitudes, rather 
than the phases. It would be interesting to extend the analysis to the non-linear scales (high Z's) for 
which the phases are non-random, and to simulations in which the initial conditions had non-random 
phases, e.g. cosmic strings and texture. An improvement to the statistic can be achieved by presenting 
the Net with much larger number of independent simulations. 

While the conventional approach of comparing rms values of the Harmonics or correlation functions 
is more straightforward and can be used to reject models, our approach which incorporates phases is 
useful when it is difficult to discriminate between models based on the amplitudes alone. As with any 
statistic, the result depends on what physical information is fed into the Net. The method could be 
generalized to present the Net with data from redshift surveys and with peculiar velocities of galaxies. 

7 Discussion 

We have presented several examples of application of Spherical Harmonics Analysis to the study of 
the large scale structure, for both cosmographical (e.g., the Puppis cluster) and cosmological inference 
(e.g., the power-spectrum and fi[]' 6 /&). We are currently developing the Harmonic expansion in 3-D 
further, correcting simultaneously (using the properties of the Harmonics) for incomplete sky-coverage, 
redshift distortion and the shot-noise by Wiener filter (Fisher, Hoffman, O.L., Lynden-Bell & Zaroubi, 
in preparation). This procedure will then be applied to new all-sky IRAS and optical redshift surveys, 
and to surveys of the peculiar velocity field. The 3-D p\ mn coefficients will allow objective non- 
parametric comparison of different surveys of light and mass in the local universe. 
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6 Observed vs. simulated Harmonics using Artificial Neural Networks 

The comparison of the observed large-scale structure with cosmological simulations remains a chal- 
lenge. The traditional methods, like the 2-point correlation function, and counts-in-cells, do not 
represent the phase information. In simple words, two realisations which have very different apparent 
features could nevertheless have the same correlation functions. Another point concerns the method- 
ology of comparing models with observations. Usually, by performing a x 2 or maximum likelihood 
test we estimate the probability for the data given a model, P(data\model). But what is really of 
interest is the probability for the model given the data, which is given by Bayes' theorem, 

P(model\data) oc P(model) P(data\model). (15) 

Artificial Neural Networks (ANN) incorporate very well the aspects of phase information and the 
Bayesian approach, without ad-hoc assumptions, e.g. that the Harmonics are drawn from a Gaussian 
density field. ANN algorithms were originally developed to model the human brain. Here we use an 
ANN model known as the Backpropagation algorithm (e.g. [7]). It consists of nodes (analogous to 
human neurons) arranged in a series of layers. In our case the input layer consists of the 24 Harmonic 
coefficients, and the output layer consists of the 3 models ('classes'). We also include 'hidden layers' 
which allow non-linearity in a complex classification space. ANNs have recently been applied to 
several problems in Astronomy, e.g. morphological classification of galaxies [26]. Non- astronomical 
applications somewhat similar to our problem are face and speech recognition and identification of 
hand-written characters. The ANN 'learns' by minimizing least-squares and fixing its free parameters, 
known as the 'weights', from a sample for which the answer is known (simulations in our case). Then 
it can predict classification for new data. The output vector can be viewed as the Bayesian a posteriori 
probability for a model given the data. Moreover, the sum of the output vector components is nearly 
unity, as expected for a probabilistic classifier. Here we train ANN on simulated Harmonics for which 
the assumed model is known, and then present the ANN with the observed Harmonics, and ask it 
to choose a favourite model, a question usually addressed in a limited way by statistics such as the 
2-point correlation function, or by visual impression. 

As an example, the observed projected (number- weighted) distribution of IRAS galaxies brighter 
than 0.7 Jy [17] was expanded in Harmonics 1 < / < 4, i.e. the galaxy distribution was compressed 
into 24 coefficients. The simulations [18] are of standard biased Cold Dark Matter (CDM) model, 
unbiased scale-free power-spectrum, k n with n = — 1 (which has more power on large scales than 
the standard CDM model) and a Poisson model. All simulations mimic the observed IRAS selection 
effects. We have used 24 simulations of each model, i.e. 72 simulations in total (i.e. assigning equal 
prior probability of a ^ for having each of the models). The Network configuration is (24;2;3), i.e. it 
includes 24 input Harmonic coefficients, 2 hidden units, and 3 output nodes representing the models. 

If presented with randomly oriented data the Net will be confused, since the individual a™'s are 
dependent on the coordinate system (only the mean-square for each / is invariant under rotation). To 
to measure relative phases, we first aligned the dipole and quadrupole of each simulation with that of 
the IRAS sample. We then carried out (O.L. & C. Scharf) the following experiments: 

(1) The ANN was trained on all 72 realisations, and was then given the same 72 realisations to 
classify. The success rate was 94%, i.e. the Net 'memorized' very well all the cases it has seen. While 
this may seem very impressive at first sight, it also indicates a possible problem. When the number 
of weights (59 in this case) is large relative to the number of data points (72 X 24 = 1728), we might 
over-fit the data. 

(2) As a more challenging task for the Net, we divided the 72 simulations into two, a training set 
of 48 simulations (16 for each model), and a testing set of 24 simulations (8 for each model). After 
training the net on the 48 simulations we presented it with the 24 simulations it had never seen before, 
the Net recognised all 8 Poisson simulations correctly, as the human eye will easily do. Of the 8 CDM 
realisations it recognized 7 correctly, and 1 wrongly as Poisson. Of the n = — 1 realisations it classified 
only 3 realisations correctly, and 5 wrongly as CDM. 



Figure 3: Reconstruction of 3-D Harmonic expansion (in redshift space) with p\ mn coefficients up to 
Harmonic l max = 15, evaluated at r = 1000 km/sec. The Local Supercluster appears here centred at 
Galactic coordinates (/ ~ 290°; b ~ 45°) and the Local Void at (/ ~ 45°; b ~ 0°). 



Bessel functions within the sphere, J Q a dr r 2 ji(k n ir) ji{k nn r) = S nn ic n l with = \a 3 {k n ia)if{k n ia), 
where c n \ is also the normalization in eq. (13). Note that the Bessel function ji(z) has an infinite 
number of zeros, which behave asymptotically like (k n ia) ~ 7r(n + //2). We therefore need to choose the 
maximal number of radial modes n max (l) for a desired spatial resolution. Other boundary conditions 
are also possible, although they are somewhat less natural, e.g. one can choose (ref. [1]) the pi m {r) 
to vanish at r = a, corresponding to the zeros of ji{k n ia) = 0. One can also choose a different set 
of radial functions (e.g. the associated Laguerre polynomials [15]), but the spherical Bessel functions 
naturally appear in our problem (see Section 4). 

The above formalism only holds for a volume-limited sample. In practice, given a flux-limited 
sample we can estimate the coefficients from the data by 

Plmn = TT~s M k nl' ^') Y* m (r t ) , 

gal ™ Tl) 

where the sum is over galaxies with r < a, and (f>(r) is the selection function. 

As an example of the method, we have applied it to the 2 Jy IRAS redshift survey [27] with 
a = 6000 km/sec, including 1888 galaxies in that sphere. The number of radial modes n max (l) was 
chosen to give a desired resolution of k n \i < 21 /a and l max = 15, resulting in 617 p\ mn coefficients. 
We show a reconstruction (in redshift space) by evaluating pi m {r) = J2 n c nlPlmnjl(k n l' r ) at r = 1000 
km/sec and reconstructing the angular Harmonic expansion up to l max = 15. Figure 3, in Galactic 
coordinates, shows the Local Supercluster and the Local Void in full glory. However, at larger distances 
the reconstruction is getting more noisy, in part because of the l/(f>(r) weighting scheme. We are 
currently exploring other weighting schemes and the incorporation of a 'Wiener filter' to improve the 
quality of the 3-D reconstruction at larger distances. 



(14) 



Figure 2: Wiener reconstruction of the 2-D 1.2 Jy IRAS galaxy sample, for Harmonics / < 15, plotted 
in Aitoff Galactic projection. The reconstruction corrects for incomplete sky coverage, as well as 
for the shot-noise. The assumed prior model is a low density CDM (with shape parameter T = 0.2 
and normalization <jg = 0.7), although a prior of standard CDM gives a very similar reconstruction. 
The reconstruction indicates that the Supergalactic Plane is connected across the Galactic Plane at 
Galactic longitude / ~ 135° and / ~ 315°. The Puppis cluster stands out at the Galactic Plane at 
/ ~ 240°. The horizontal dashed lines at b = ±5° mark the major 'Zone of Avoidance' in the IRAS 
sample. 

Preliminary application to simulated samples with larger incompleteness indicate bias in the 
method, which we are currently exploring. The method can be extended to 3-D in both real and 
redshift spaces, and applied to other cosmic phenomena such as the COBE Microwave Background 
maps. 

5 3-D orthogonal Harmonic expansion 

An obvious extension of the above methods is to three dimensions. We expand the fluctuations in the 
density field in Spherical Harmonics Y\ m and Spherical Bessel functions ji(z) (cf. ref. [1]): 

P( r ) ~ P = ^2^2 Plm{r)Yi m {r) = J2 Cnl P lmn ji( k nl' r ) Y lm (r) , (13) 

l m I m n 

for / > 0. Similarly, the fluctuation in the potential can be expanded as $ = J2iJ2 m ^lm( r )Yi m (r), 
and the two quantities are related by Poisson equation, V 2 [$/ m (r)Y/ m (r)] = — 4irGpi m (r)Yi m (r). Let 
us assume that the data are given within a sphere of radius a, such that inside the sphere the desired 
density fluctuation is specified by pi m (r), but for r > a the fluctuation is pi m {r) = (this simply 
reflects our ignorance about the density field out there; the fluctuations do not vanish of course 
at large distances). Hence, inside the sphere pi m {r) cx k 2 ^i m (r) cx ji(kr), and outside the sphere 
^l m (r) oc r~(' +1 ). It is sensible to require the potential and its logarithmic derivative to be continuous 
at r = a (D. Lynden-Bell, private communication). It then follows that the discrete A;'s are the zeros of 
jl-i(k n iia) = 0, where /' = /—!. It can be shown that this condition ensures the orthogonality of the 



4 Wiener filter for Spherical Harmonic reconstruction 



The analysis of whole-sky galaxy surveys commonly suffers from the problems of shot-noise (due to 
the discreteness of objects) and incomplete sky coverage (e.g., at the Zone of Avoidance). Here we 
discuss a method [14] for correcting for these effects using a Bayesian framework, and the orthogonality 
property of the Harmonics. The recovery of complete distribution from noisy and incomplete data is a 
classic problem of inversion. It is well known that a straightforward inversion is unstable, and hence a 
regularisation of some sort is essential. In the Bayesian spirit we are using here raw data and a 'prior 
model' to produce 'improved data'. 

We formulate our problem as follows: What are the full-sky noise-free Harmonics a\ m given the 
observed Harmonics c; mj0 6 s , the mask W, and a prior model for the power-spectrum of fluctuations? 
The observed Harmonics (with the masked regions filled in uniformly according to the mean) are 
related to the underlying 'true' whole-sky Harmonics by (cf. [20], eq. 46.33) 

Clm,obs = J2J2 W W m 'W lW + (Ta] (9) 

where the monopole term (/' = 0) is excluded, and we have added the shot-noise for number-weighted 
ai m Harmonics with variance (a^) = Af (the mean number of galaxies per steradian, independent of / 
in this case). 

By the rule of conditional probability and the assumption that the density fluctuations are drawn 
from a Gaussian random field we can write 

Pf n \„ \ -Pg( a 7 C obs) nn s 

P{a\c obs ) = — — — — , (10) 

rG(Cobs) 

where the vector c 6 s represents the set of observed Harmonics {c; mj0 6 s } and Pq stands for a Gaussian 
distribution function with variance and covariance which depend on an assumed power-spectrum. This 
is a special case of constrained realizations formalism ([8] and Y. Hoffman in this volume), but here 
the formulation and computation are greatly simplified due to the orthogonality of the Harmonics. 
The full derivation of maximizing the probability (eq. 10) with respect to a will be given elsewhere 
(Zaroubi et al. , in preparation), but the answer for the 'mean field' reconstruction is simply 

a = FW 1 c obs , (11) 

with 

F = dtag{—^—}. (12) 

One can also express the scatter around this value. It can be shown that this result also gives the 
minimum of the variance (|a — a| 2 ) for a desired filter matrix F. This F matrix is in fact the well-known 
Wiener filter (the ratio of signal to signal+noise) commonly used in signal processing (e.g. [22]). Note 
that it requires a priori knowledge of the variances in both the signal and the noise. When the noise 
is negligible the factor is approaching unity, but when it is significant the measurement is attenuated 
accordingly. Even if the sky coverage is 4ir (W = I), the Wiener filter can be used to remove shot- 
noise, giving the most probable picture of the underlying 'continuous' density field. In the case of full 
sky coverage, only the amplitudes are affected by the correction, not the relative phases. For example, 
the dipole direction is not affected by the shot-noise, only its amplitude. Of course, if the sky coverage 
is incomplete, both the amplitudes and the phases are corrected. The reconstruction also depends on 
how many Harmonics are observed (l max ) and how many are desired for the reconstruction (l' max ). 
Note also that the method is non-iterative. 

Here we apply the method to the 2-D sample of 5313 IRAS galaxies brighter than 1.2 Jy ([4], 
[5] ), covering 88 % of the sky. Practically, rather than inverting W (eq. 11) we solve the equation 
c t, s = WF _1 a using the Singular Value Decomposition algorithm. The 'best reconstruction' is shown 
in Figure 2, indicating the connectivity of the Supergalactic Plane where crossed by the Galactic Plane, 
and confirming the Puppis cluster. 
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where the 'window function' which depends on the radial selection function (f>(r) is 

*P(k) = [dr r 2 <j>{r)f{r)jikr) . 



(4) 



We have assumed that the power-spectrum of the galaxies is b 2 P m (k) ('linear biasing'). To this 
'cosmic scatter' one should also add the 'shot noise' due to the discreteness of objects, {\ai m \ 2 ) sn = 
Jdrr 2 ^r)f 2 (r) 

3.2 Galaxy Harmonics in redshift space 

The power-spectrum derived from a redshift survey will differ from the one in real-space (Kaiser[9]). 
This is also seen in the Harmonic analysis of redshift surveys (Scharf & Lahav[25]), and the Harmonics 
in redshift space can be formulated (Fisher, Scharf & Lahav[6]) for a flux- limited survey in linear theory 
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where the redshift distortion window function is 
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For / = 1 (dipole) this recovers Kaiser's 'rocket effect' [9]. In the special case that both (f>(r) and f(r) 
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are power laws, the real and redshift Harmonics are simply related by 
where T\ depends on the power-indices of the selection and weighting functions, and (3 = J7q' 6 /6. 

We have applied (ref. [6] and K. Fisher in this volume) eq. (5) and a Maximum Likelihood 
analysis to the 1.2 Jy redshift survey ([4]). If we fix the normalization of IRAS galaxies in real space 
(og = 0.69 ± 0.04) and solve for the shape parameter (T = Clh) of a CDM power-spectrum, we get 
r = 0.17 ± 0.05 and (3 = 0.94 ± 0.17 (1 — a). This value of the shape-parameter is in accord with other 
measurements of excess power on large scale (e.g. [5]), and (3 for IRAS (which is strongly coupled 
with eg ) is again found to be higher than (3 deduced from optical studies (e.g., [16], [12], and M. 
Hudson in this volume). It is interesting to note that although we have not used here any independent 
measurement of peculiar velocities (i.e. distance indicators), our result for (3 is in agreement, and of 
comparable accuracy, to the results deduced by comparing the peculiar velocity field (POTENT) with 
IRAS redshift surveys (e.g. [2]). 

3.3 Peculiar Velocity Harmonics 

The Harmonics for the peculiar velocities are simply the difference between Harmonics in redshift 
space and real space: 
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but in the case of incomplete sky coverage the window function is more complicated [21]. 
3.4 The Microwave Background (Sachs- Wolfe) 

Finally, the much-discussed Sachs- Wolfe fluctuations in the Microwave-Background are written in 
terms of Harmonics as: 
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Figure 1: Spherical Harmonic reconstruction with coefficients up to l max = 10, of galaxies with 500 
km/sec < czlq < 3000 km/sec in the 2Jy IRAS (\b\ > 5°) combined with our Puppis sample (\b\ < 5°). 
Plots are equal area hemispheres with the left-hand side plot centred on Galactic / = 240°, b = 0°. 
The Galactic Plane runs horizontally across the plots (solid line), dashed lines bound the region 
\b\ < 5° and longitudes are indicated. The South Galactic hemisphere is at the top. The lightest 
solid contour is at the mean, dashed and solid contours indicate densities below and above the mean 
respectively. Contour separation is 3 times the shot-noise level. Associations with local structures 
are labelled: V - Virgo, PU - Puppis, F - Fornax E - Eridanus, Leo - Leo, Cen - Centaurus, PIT - 
Pavo-Indus-Telescopium, Cam - Camelopardalis and UM - Ursa Major. From ref. [13]. 

studied independently by other groups ([10] and [28]). To quantify the importance of Puppis relative 
to other structures in the Local Universe we supplemented the IRAS 2 Jy \b\ > 5° redshift survey [27] 
with an IRAS-selected sample in the direction of Puppis (\b\ < 5°; 230° < / < 260°), which consists 
of 32 identified galaxies, 12 of them with measured redshift. We find [13] that the projected number 
counts of galaxies brighter than 2 Jy in Puppis is about half that of Virgo. A Spherical Harmonic 
reconstruction (Figure 1) shows that out to a distance of 3000 km/sec Puppis is second only to 
Virgo. We estimate that Puppis (which lies below the Supergalactic Plane) may contribute at least 
30 km/sec to the motion of the Local Group perpendicular to the Supergalactic Plane (Vsgz = —370 
km/sec in the Cosmic Microwave Background frame). Together with the Local Void (above the 
Supergalactic Plane) and Fornax and Eridanus (below it) this may explain the origin of the so-called 
"Local Anomaly". 

3 Spherical Harmonics as probes of the power-spectrum 

The mean-square of Harmonics can be related to the power-spectrum of mass fluctuations P m (k) = 
{\Sk\ 2 ) in Fourier space for a variety of cosmic phenomena. A useful identity in deriving the relations 
below is the expansion of a plane wave in spherical waves, e* k ' r = 4ir J2i Em (0' jl(^ r ) Yi m (r) Y/_ m (k) , 
where j\ is the spherical Bessel function. 

3.1 Galaxy Harmonics in real space 

If the distances to galaxies are known then an estimator for the Harmonic coefficients is a^ m = 
J2gal /( r «)^ / /m(^i)' wnere f( r ) ' ls a 'weighting function'. The mean-square prediction is ([24], [25]): 



where the radial coefficients a/ m (r) are estimated using an appropriate choice of weighting function 
f(r) and the orthogonality property of the Harmonics by 

aim(r) = '£f(r i )Y l * m (r i ), (2) 

gal 

where the sum is over the observed galaxies. Spherical Harmonic Analysis (SHA) has been discussed 
for analysing projected surveys about 20 years ago [19] but was not that useful given the poor sky 
coverage of the samples existing at the time. More recently SHA has been reconsidered and applied 
to IRAS surveys ([3], [24], [25], [13], [6], [23] ) and to the peculiar velocity field ([21], [15], [11]). 
Some of the advantages of SHA in probing the large scale structure are: 

(i) The Harmonics retain both amplitude and phase information, hence providing a unified language 
for both cosmography and statistics. 

(ii) A spherical orthogonal coordinate system is the natural one for analysing whole-sky distribu- 
tions. 

(iii) The SHA provides natural smoothing with angular resolution ~ it /I. It covers a wide range 
of scales, from the dipole to the small scales usually probed by correlation functions. 

(iv) The SHA is an efficient data compression procedure. The compressed data can then be used 
as input to various statistical tools such as Maximum Likelihood, reconstruction methods, or Artificial 
Neural Networks (see below). 

(v) The estimated coefficients can be used for comparison of different tracers of the density field 
(e.g. IRAS and optical), as well as with the peculiar velocity velocity field and Background radiations. 

(vi) The much-discussed gravity dipole (/ = 1) and quadrupole (/ = 2) are only special cases of 
the SHA. 

(vii) At small angles the SHA becomes the familiar 2-D Fourier transform. The mean-square of 
amplitudes is simply related to the 3-D power-spectrum or the 2-point correlation function. 

One might argue that SHA assumes that the observer is at the centre of the universe and is 
therefore a biased statistic. However, shot-noise increases with radial distance from us in flux limited 
surveys. Moreover, the analysis can be done from other centres as well. Another major problem of the 
SHA is incomplete sky coverage, in particular due to obscuration and/or confusion by the Galactic 
Plane. This can be accounted for by modelling an angular window function([19], [24]), or by a Wiener 
reconstruction ([14] and below). 

Here we discuss several examples which illustrate the utility of the Harmonics in overcoming some 
of the key problems in the study of the large scale structure: shot-noise, redshift distortion and 
incomplete sky coverage, and in estimating properties such as the power-spectrum, biasing and the 
density parameter fio- We summarize studies on a projected IRAS sample of galaxies brighter than 0.7 
Jy [24], the 2Jy IRAS redshift survey [25], and the Puppis cluster of galaxies behind the Galactic Plane 
[13]. We then review the Harmonics as probes of the power-spectrum and cosmological parameters, 
in particular using the SHA to analyse redshift distortion [6] and to estimate the combination of the 
density and bias parameters J7q' 6 /6. Then we discuss on-going projects of Wiener filter reconstruction 
of noisy surveys with incomplete sky coverage [14], a 3-D orthogonal expansion and pattern recognition 
with Harmonics and Artificial Neural Networks. Further details and discussion on Harmonics are given 
in this volume by K. Fisher (on redshift distortion) and by C. Scharf (on cross-correlation of optical 
and IRAS Harmonics). 

2 Cosmography with Spherical Harmonics: the Puppis cluster 

The expansion in Spherical Harmonics can be used to reconstruct the surface density or brightness 
up to a certain Harmonic l max - Various reconstructions by Scharf et al. ([24], [25] [13] and [23]), 
reveal familiar and new structures, and illustrate for example the 'tug of war' between the Great 
Attractor/Centaurus and Perseus-Pisces superclusters. Figure 1 in ref. [24] indicated a significant 
overdensity in the direction of Puppis (/ ~ 240°; b ~ 0°). This cluster has also been noticed and 
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Abstract 



The orthonormal set of Spherical Harmonics (Y; m ) provides a natural way of expanding whole- 
sky redshift and peculiar velocity surveys. This decomposition retains both amplitude and phase 
information, hence providing a unified language to describe the local cosmography as well as the 
power-spectrum of the galaxy density and velocity fields. We apply the method to IRAS projected 
and redshift surveys, to investigate redshift distortion and the value of Slo, and to explore the Puppis 
cluster hidden behind the Galactic Plane. We further discuss the choice of radial function and the 
estimation of the coefficients in 3-D, and the removal of shot-noise and correction for incomplete 
sky-coverage by a 'Wiener filter'. The data compressed by the Harmonic decomposition can then 
be analysed by Artificial Neural Networks in comparison with simulations. 

1 Introduction - why Spherical Harmonics ? 

The rapid progress made in observations of large scale structure and in particular the production 
of 'nearly all-sky' redshift surveys, calls for new approaches to quantify the galaxy distribution and 
peculiar motions. The much-discussed statistics such as correlation functions and 'counts-in-cells' are 
most useful for estimating the underlying power-spectrum of the density field, but ignore the phase 
information, which is so crucial for describing features such as the Supergalactic Plane. 

In the light of new IRAS and optical 'whole-sky' redshift surveys and the intensive work on pe- 
culiar velocity surveys it is natural to discuss the density and velocity fields by means of Spher- 
ical Harmonics. In brief, the density field is expanded using the orthogonal set of Harmonics, 
Yi m (Q,<f>) oc P ; ' m ' (cos#) exp(im^)), where 6 and <f> are the spherical polar angles, and P ; ' m ''s are 
the associated Legendre Polynomials of degree / and order m. Generally a 3-D density field can be 
expanded out to l max by 





m = — l 



